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ABSTRACT 


A general computer program has been developed to 
optimize an n-stage vapour-compression refrigeration system 
based on minimum total cost comprising initial and running 
expenditures . The condensing temperature has also been treated 
as a decision variable and the costs of all the compressors and 
the condenser and the cost of cooling water have been included 
in the total cost. Optimum interstage pressure and condensing 
temperature have been determined for two-stage R-12 and R-22 
systems, making use of available compressor and condenser sizes 
and their costs. Effects of subcooling, superheating, approach 
and compressor efficiency have been considered to make the 
results more relevant. Results have been presented for various 
design ambient temperatures and also considering year-round 
variations in ambient temperature. Optimum interstage pressures 
obtained by the present method have been compared with data 
available in the literature and has been found to be greater 
than that of all other investigators. It has been found to 
be about 15% higher than the geometric mean pressure. 



CHAPTER - 1 


INTRODUCTION 


1.1 DESCRIPTION: 

It is well known that the Carnot refrigeration system 
having two reversible adiabatic and two reversible isothermal 
processes is the perfect system of refrigeration. Its coeff- 
icient of performance (COP^) depends entirely on the two 

temperature limits. But/ in practice it is not possible to 

( 1 ) 

operate a Carnot cycle. Yet it is desirable that any 

actual refrigeration system should have its coefficient of 
performance (COP) close to COP^. It is sought to be achieved 
with the help of a phase change cycle consisting of an iso- 
thermal evaporation/ a reversible adiabatic compression/ an 
isobaric condensation having a very small non-isothermal 
process, and an isenthalpic expansion. The COP of this system 
tends to COP^, if the difference between the condenser and the 
evaporator temperatures is small. But, in practice the con- 
densing temperature is decided by the ambient temperature ar^ 
the evaporator temperature is fixed by the cooling require- 
ments. For low evaporator temperature, the single-stage 
vapour compression system referred to above has many disadvan- 
tages: 

(i) Increase in compression ratio reduces the volumetric 
efficiency of the compressor - 1 - 

- 1 It means that larger compressor 
displaceirent would be required for a specified capacity 
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(ii) Leakage past the piston is increased due to large 

pressure difference. Also there is overheating of the 
compressor body giving rise to lubrication problems, 
and the overall efficiency of the compressor is 
considerably reduced. 

(iii) Due to high pressure ratio across the expansion 

valve more flash vapour is produced after the expansion 

which reduces the refrigerating effect- On the other 

hand the work of compression is increased, reducing COP. 

(iv) Smaller the COP, larger is the heat rejection in the 

condenser [Q = Q ,• (1 + 1/COP )] . This 

'-rejected cooling ' -* 

results in increased size of the condenser causing 
higher cost . 

To overcome these difficulties, multistage (generally 
2- or 3-stage) vapour compression system is generally used. In 
such a system the vapour from the evaporator is compressed to 
the condenser pressure in two or more stages with intercooling 
in-between. An n-stage system is exhibited in Figure 2.1 along 
with p-h diagram. A 2-stage system is shown in Figure 2.4. 

As compared to a single-stage system, the additional components 
in a 2-stage system are: a compressor, an expansion valve, 
and a flash intercooler. Liquid refrigerant together with 
some flashed vapour enters the evaporator (E), absorbs heat 
and gets converted into vapour. The vapour is then compressed 
by the LP compressor to an intermediate pressure (and not to 
the condensing pressure as is done in a single-stage system), 
at which the flash intercooler (FI) is maintained. The comp- 
ressed vapour gets desuperheated by the liquid refrigerant in 
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the flash intercooler. The heat released due to desuper- 
heating augments the vapour flow rate through the HP compre- 
ssor. This vapour, the flash vapour due to throttling of the 
condensate and the LP saturated vapour are then compressed to 
the condenser pressure by the HP compressor where heat transfer 
to the cooling medium causes condensation of the vapour. The 
condensate is throttled upto flash intercooler pressure through 
the expansion valve (EVl) . The saturated liquid from the flash 
intercooler is supplied to the evaporator through the expansion 
valve (EV2)/ completing the cycle. 

Referring to Figures 2.3 and 2.5, it can be observed 
that the throttling process is closer to the saturated liquid 
line. It implies therefore that the refrigeration effect per 
unit mass is increased . On the other hand , on the compression 
side, the processes are closer to the saturated vapour line. 

It renders saving in work of compression and also less super- 
heated vapour which means less heat rejection. Thus the 
limitations of a single stage-system are adequately overcome. 
Further, a multistage system also has some additional advan- 
tages: 

(i) By proper arrangement of the cranks for all the stages, 
torque can be made more even; thus requiring a lighter 
flywheel. 

(ii) Here LP compressor is subjected to lower pressure and 
so the cylinder wall thickness is smaller. Thus a 
saving in material cost occurs. Also, only the HP 
compressor needs to be made of special materials to 
withstand the high pressuure. 
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But, in this system the number of auxiliaries and 
spare parts is more, giving rise to increased initial invest- 
ment. However, this disadvantage is offset by reduced running 
costs for higher capacities and low temperature applications. 

1.2 LITERATURE REVIEW: 

The performance of a multistage system is dependent 
upon the interstage pressure (s) or temperature ( s ) selected. 
There exists a set of such pressure (s) which gives (give), the 
best perfomance of the system and it (they) is (are) called 
"Optimum interstage pressure ( s )" . It (they) is (are) obtained 
by maximizing or minimizing certain performance criterion 
(objective function) and this criterion must be explicitly 
spelt out. So far three optimization criteria have been 
envisaged. They are: 

(i) minimization of work input per Kg of refrigerant 
through .the evaporator^ ^ 

( 8 ) 

(ii) maximization of the COP of the system^ ^ 

(iii) minimization of power per ton of cooling. 

The last two criteria are basically equivalent to 
each other and so give the same results. But the work input 
per Kg of refrigerant f low is not a complete indicator of the 
performance of the system because it does not include the 
variation in cooling produced. 

In general for n- stage compression of perfect gas 
with perfect intercooling, the' interstage pressures are in 
geometric progression given by 

= P^o/P,-i = . . - pVp,-«._i = (P.,/Pi)^^" 
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for minimum work of compression and, in particular, for a two- 
stage compression the interstage pressure (p^) is given by 

Pi “ (PhPi>°’" 

where p^ is the discharge pressure of the HP compressor and 
Pj^ is the suction pressure of the LP compressor. But p^ given 
above (eqn. (1.1b)) has been shown to be lower than that for 
maximum COP . ^ ^ ^ 

( 2 ) 

Soumerai has developed an iterative scheme which 

uses a graphical method to find the interstage pressure. 

. ( 3 ) 

Keshwani and Rastogx obtained the optimum interstage pres- 
sure corresponding to minimum work. They assumed the mass of 
the flashed vapour to be invariant with the intermediate 
pressure. This assumption is unrealistic. Also the pressure 
at which the mass of the flashed vapour evolves is not recomm- 
ended. Surprisingly their results are found to be unrealistic 
for R-12 and R-22 systems. 

( 3 ) 

The work of Keshwani and Rastogi has been extended 
( 4 ) 

by Ram Lai . ^ He accounted for the variation in mass of 

flashed vapour with intermediate pressure. But this method 

involves the generation of complete superheated properties 

which is (quite cumbersome. Also it corresponds to minimum 

(5) 

work. Verma and Charan • used linear variation in mass of 
flashed vapour with interstage pressure, rendering simplifici- 
cation in Ram Lai’s work. Their expressions consist of a set 
of differential equations. But unfortunately any numerical 
result based on their analysis to check the method was not 
established. Their method is also based on minimum work. 
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Surprisingly the interstage pressure reported corresponds to 
maximum COP . 

Gosney^^^ has studied two-stage vapour — compression 

systems, but only for a few operating conditions. Arora and 
( 1 ) 

Dhar have analysed n-stage system and reported results 
based on the minimum work criterion. Moreover, they have 
considered the ideal cycle without the effects of compression 
efficiency, subcooling and superheating. Such a cycle is an 
unrealistic representation of the actual system. Further, 
one has to resort to computer for calculation of optimum inter- 
stage pressure for each set of operating conditions. 

/ Q \ 

Prasad^ has developed a correlation for a two-stage 
R-12 system in terms of the condenser and evaporator temper- 
atures. COP has been considered as the objective function 
and an approximate method, neglecting the superheat horn, has 
been used. It is very well valid for R-12 system because it 
has a very small area of superheat horn as compared to work 
of compression. The method completely avoids the use of 
superheated properties which is a great simplification without 
significant loss of accuracy. But it cannot, unfortunately, 
be applied to other systems because of a significant area of 
superheat horn as compared to the work of compression. 

1.3 PRESENT STUDY: 

In the past investigations , attention has been focussed 
on the thermodynamic performance of the system. In the present 
study the focus has been shifted to cost aspects of the system 


as a whole 
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Wnen we change from a single-stage system to a 
multi-stage system^ the number of compressors, intercoolers 
and expansion valves increases . Optimization of the thermo- 
dynamic performance overlooks the costs of these components, 
and hence the calculation of optimum interstage pressure is 
more of an academic interest rather than that of commercial 
requirements. To make the result more relevant, the costs of 
these components has also been included. Electricity is 
consumed in running the compressors and cooling water is 
required for the condenser. So the costs of .these items 
have been considered besides the costs of both the compressors. 

Evidently for maximum COP the condensing temperature 
should be as low as possible and due to natural limitations 
it may be kept as close to the ambient temperature as possible. 
But as the condensing temperature is reduced, the mass flow 
rate of cooling water increases, giving rise to increased 
cost. For a two-stage system, we have thus both the inter- 
mediate temperature and the condensing temperature as decision 
variables leading to a 2-dimensional minimization problem. In 
the present work, rather than fixing the condensing temperature, 
an optimum condensing temperature has been obtained along with 
optimum cooling water rate and optimum interstage pressure 
with respect to minimum cost. 

In the present study, effects of subcooling after 
condensation and superheating in the evaporator have been 
considered. Also the effect of variation in approach temper- 
ature (difference between condensing and outlet cooling water 
temperatures) on the optimum quantities has been considered. 
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The effects of tonnage capacity of the system has also been 
considered. Optimum quantities have been obtained for various 
evaporator temperatures. Also all the optimum quantities have 
been obtained for different design ambient conditions. Thus, 
for a known set of design parameters (viz. subcooling, super- 
heating, approach, the evaporator temperature and the tonnage'' 
and for any given design ambient temperature of a given loca- 
tion, an optimum system can be recommended which includes size 
of the compressors for HP and LP side condenser size and 
economic cooling water rate . 

As the compressors and other equipments are available 
only in certain discrete sizes, the cost function comes out 
to be a discontinuous function of the decision variables. A 
general computer program has been developed to minimize a 
function of n variables. Optimization has been done for two- 
stage R-12 and R-22 systems as the cost data for other refri- 
gerant compressors were not received from the sources. 

The work has been extended to include the effect of 
variations in ambient temperatures. In the absence of any 
control system, the condensing and the inter-stage temperature 
would change as the ambient temperature changes. Hence the 
system would deviate from optimality. But the prediction of 
these variations is not easy. In the present work, rather 
than optimizing for any one ambient temperature, a joint ef £€>-"■ 
of monthly mean maximum temperatures for a particular place 
has been optimized . Kanpur has been selected as the place 
and the temperature data have been taken from reference (10). 
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For the present investigation, the costs of various 
components of the system have been obtained from Messrs Frick 
India Ltd., New Delhi, Messrs United Refrigeration Works, 

Kanpur and Messrs Kaiyan Cooling Corporation, Kanpur. Cost of 
water has been obtained from Jal Sansthan, Kanpur. Cost of 
electricity was obtained for various years from KESA House, 
Kanpur arxi a logistic curve has been fitted through this past 
data to forecast costs for the future. Properties of refri- 

( 13 '■ 

gerants have been used from correlations obtained by Prasad. ' 
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CHAPTER - 2 


PROBLEM FORMULATION 


2.1 MAIHEMATICAL MODELLING: 

2.1.1 n-stage System: 


Figure 2.1 shows a generalized system consisting of 

n-stage compression with flash intercooling. In this system 

the evaporator vapour is compressed through successive stages 

with desuperheating in the flash intercoolers until the vapour 

is discharged into the condenser. To derive generalized results, 
tin 

the i stage as shown in Figure 2.2 is considered. xMass and 
energy conservation for the flash intercooler of this stage 
yields : 


■^i ~ ^i~1^^2i-2 ” ^4n-2i+3^'^^^2i " ^4n-2i+l^ 

for i = 2, 3/ .../ n 


( 2 . 1 ) 


* til 

where m^^ is the mass flow rate through the i -stage compressor. 

Also, 


A 


, = 3.5 0^/(h, - 


( 2 . 2 ) 


. th 

Metric Horse Power required by the i compressor is: 


‘‘h = - h2^.i)/0.736 


(2.3) 


, th 


If ^ be the compression efficiency of the i stage, one can 

get y: ; ^ • : 


^2i “ ^2i-l ^^2i‘ " ^2i-l^^ 


Cl 


(2.4) 
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where is the enthalpy at the end of isentropic compression. 

For the entire system, 

n 

COP = 4.75543 6 / 5 " HP, (2.5) 

c i^l 

Heat rejected in the condenser is : 


"'h 


III (h.., 
n 2n 


" ^2n+l^ 


( 2 . 6 ) 


Also, 


n 


Q,. 


= 3.5 Q + 0.736 21' 


i=l 


(2.7) 


If be the temperature rise of the cooling water in the 

condenser, the mass flow rate of cooling water is: 


ro = Q. /(c -AT) 
w ^ pw ^ w^ 


( 2 . 8 ) 


needed in this expression is calculated on certain assump- 
tions as explained below: 

Referring to the real condensation process shown in 
the T-s diagram in Figure 2.5, we observe that the condensation 
process consists of three parts: 

(i) desuperheating of vapour, process A-B (in actual 

practice it occurs in two parts — (a) desuperheating 
and (b) combined desuperheating and corxiensation) 

(ii) condensation of vapour to liquid, process B-C 

(iii) subcooling of the condensate, process C-D. 

The total heat rejection is given by the area A-B-C-D-d-a-A, 
consisting of heat rejections in the above three processes. But 
it is extremely difficult to predict the zone near point B where 
both condensation and desuperheating continue having extremely 
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complicated heat transfer process. On the other hand, if we 
consider the entire condensation to be occurring at a fixed 
temperature the heat rejected will be given by the area 
E-F-d-a-S which is smaller than the original area by an amount 
given by the area ASB minus the area CFD. But this difference 
is not appreciable in comparison to the total area and if 
required it may be neglected in the light of a major simplifi- 
cation. Validity of this approximation is also supported by 

(15) 

other researchers. However', the computation has been 

carried out for the same heat rejection. Only the condensing 
temperature has not been adjusted due to expected negligible 
change in So with this approximation, referring to Figure 

2 . 6 , we see that the exit water temperature is given by 

T = T - T 
WO h ap 

and so. 




T - T . 

WO W1 


T. - T . - T 
h wi ap 


(2.9) 


2.1.2 Two-stage System: 

For a two-stage system (Figure 2.4), we have n = 2 and 
i = 1 and 2, For these values of n and i one gets from eqns. 


(2.1)-(2.7), 

= 3,5 “ ^7^ (2.10c 

^2 = ^1(^2 “ ^7^^ 

HPl = m^(h 2 ^ 

HP2 = 1112(^4 - h2)/0.736 (2.10d 
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+ (h^, - (2.1P-e) 

= h3 + (h^, - h^)/r^^.^ (2.icf; 

COP = 4.75543 6^/(HPl + HP2 ) (2.10^j 

= 3.5 + 0.736(HP1 + HP2) {2. irk) 


is given by eqn. ( 2 . 8 ) itself. 

Correlations are available for saturated properties and specific 
heats (for finding enthalpies of superheated states) as a 
function of saturation temperature.^ Thus, we can find all 

the enthalpies ; 

h . = h , , + c T . T ( 2 . 1 j. ; 

1 1 ' pi s 

s^ = s^, + Cp^ ln[l + Tg/(T^ + 273.15)] (2.12) 

h^, = h 3 + Cp 2 (T^ + 273.l5)[exp| (s^ - S 3 )/Cp 2 ] - l] (2.13- 

^ 4 < = ^ 4 " %4^’^h 273.15)[exp I (S 3 - O/c^^j- - l] (2.1^^’ 
^5 " '^c'^ (2.15) 

where, 

= evaporator temperature 
= condenser temperature 
T = superheating in the evaporator 
= subcooling of the condensate 
Tj^ = interstage temperature . 

and i^own functions of pressure ratios P 2 /l^ 

P 4 /P 3 respectively and have been taken from reference ( 12 ) as ; 

1/n^ = 1.13623 + 0.113289 r - 0.0334529 r^ 

+ 0.486757 X lo"^ r^ - 0.2134 x lO"^ r^ (2.16; 
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where r is the pressure ratio. The pi'essures are themselves 

{ 13 ) 

tunctions ot saturation temperatures. Thus all the above 

expressions given in eqns. (2.1)-(2.16) can be expressed in 
terms of temperatures. An explicit expression has not been 
attempted at as it becomes extremely involved. Moreover ^ it 
serves no useful purpose because we can directly use tne above 
implicit expressions in computer program to evaluate the 
objective function. 

Thus, we finally find the following functional rela- 
tionships : 


HPl 

= HP14T.. T^, , T^, 6^4 

(2.17c. 

HP 2 

= Th-- '^wl' Tap- T 3 , T^, 6^4 

(2.17b'; 

Qh 

= T„., Q^4 

(2.170) 

m 

w 

= VA- Tj,; T^- V- T^p, T^, Tp, ip 

( 2.17c. • 

COP 

= C0P4T^, T^, T^, ip 

( 2 .l 7 e' 

2.2 

OBJECTIVE FUNCTION: 


2 . 2 . 

1 For Design Ambient Conditions 



In the present work the objective function considered 
is the average annual cost for the entire life span in teirms 
of the present value . Thus the total cost incurred (in terms 
of the present value) over the entire life span of L years ha'' 
been divided by L to compute the objective function. It 
consists of initial and running costs. 
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During any particular optimization, the evaporator 
size is fixed. So only the costs of all the compressors and 
the condenser has been included in the list of initial costs . 
The running costs include the cost of power (electricity), the 
cost of cooling water and the cost of maintenance. 

The costs of the compressor and the condenser are 
functions of HP required for the compressor and heat rejection 
(Qj^) in the condenser respectively. Thus, 

n . . 


^initial = L ^comp'^^V ^cond-^V3 


(2.18) 


For the cost of electricity a logistic curve of the 
form a/[b + exp(-cx)] has been fitted through the past data 
and future costs have been forecast. Thus the cost rates of 


iSt ^nd 

c . in the 1 ,2 

eL 


electricity (Rs/kWh) b c _, ..., 

0 X 0 X 

years can be obtained. Also the cost rate of the cooling 

3 

water c (Rs/m ) has been obtained from "Jal Sansthan, Kanpur" 
w 

If R be the rate of interest, the present values of 
these costs can be computed as; 


'e ,eff active 


1 — 1 i-1 

i > c ./(I + R)-^ ^ 
L er 


(2.19a) 


1 i-1 

= r c,/(i + R)"" 

^ 1=1 « 


w, effective 
So , 

C = 6447 .36 * OF c 

power 


n 

* 5 . 


e, effective i 


(2.19b) 


( 2 . 20 ) 


and 


"water 


31536 * OF * c 


w,eff active 


m 


w 


( 2 . 21 ) 
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in Rs/year, where OF is the part of the day for which the 
system has to operate and is called operating factor. 

Maintenance cost has been considered as 10% of the 
equipment cost based on data given in reference (14). Thus, 


maint 


= 0.10 * 


'^initial 


L 


[(1 + R)" - 1]/[R(1 + R)^'"] 


( 2 . 22 ) 


So j 


c — c. 4-r -i- c ('? 7 

running power water maint . 

The objective function (total annual cost) is then given by: 


Z = 4- C 

■"initial running 


(2.24) 


and its functional dependence is expressed by 


Z = Z4T., Tj,; T„i, , T^, Q^, , C^, 

(2.25) 

Our aim will be to minimize this Z subject to the constraints: 


Tj_ ^ \ ( 2 . 26 ) 

and , 



These relations express the fact that the interstage temper- 
ature must lie between the evaporator .and the condensing 
temperatures and that the condensing temperature must exceed 
the inlet cooling water temperature plus the approach . 

2.2,2 Variations in Ambient Temperatures; 

In the previous section, the objective function has 
been formulated considering the ambient temperature as a 
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design parameter. Thus the cost can be minimized for any- 
selected design ambient temperature. When there is no 
control mechanism to control the condensing and interstage 
temperatures at their optimal levels then the effect of vari- 
ations in the ambient temperature should be incorporated in 
the analysis . 

As the ambient temperature varies throughout the year, 
it is not easy to predict the deviations of and from 
optimal levels . In the present study , rather than optimising 
the objective function for any single ambient temperature, a 
joint effect of the monthly temperature variations has been 
optimized . It has been still assumed that and are not 

varied -when the ambient temperature changes. Thus if 

'^wl2 monthly average -water temperatures for 

the months of January, February, ..., December, then the mass 
flow rates of cooling water corresponding to these temper- 
a-cures will be 


m 


w j 


- Aj - 


(2.28) 


We consider an average of all these water flow rates and take 


Q 


m 


h 


w- , av 


12 

s;. (: 


12c ^T - T - T 

^"^^pw j=l 'h wj ap 


(2.29) 


We use this m to compute the cost of water. This proce- 
w,av 

dure amounts to a kind of randomization over the ambient 
temperatures. Now we have 


'water 


31536 * OF * c 


'w,eff active 


m 


w ,av 


(2.30) 



is 


and C. ...^wC , C . ^ are still given by the expre- 

rnitial' power maint. j 

ssions (2.18)/ (2.20)/ (2.22), respectively. Also the final 
objective function Z is given by the same expression (2.24) 


and its functional dependence now becomes: 


Z - Z4T. , T^; T^, , c^, 


(2.31) 


Thus this new modification only changes the cost of 
water. With this approach we do not get an optimum design 
for any particular ambient conditions/ but the total cost for 
the year as a whole will be minimized. 
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CHAP'TER - 3 

OPTIMIZATION 


3.1 NATURE OF THE FuICTION: 

Our aim is to minimize the cost of the system expre- 
ssed by eqns . (2.25) and (2.31). The first step in a syste- 
matic approach to the solution of the problem requires a 
knowledge of the nature of the function and accordingly the 
technique to be adopted is to be selected . 

The objective function Z is a function of two vari- 
ables T^ and T^ for a given set. Variations in Z with either' 
of the variables / keeping the other variable fixed, are shown 
in Figures 3.1a and 3.1b for typical cases. The jumps appe- 
aring in these plots correspond to a change in one of the 
compressor capacities as compressors are available only in 
certain discrete sizes. Whenever a bigger or a smaller size 
compressor is required, there is an increase or decrease in 
the initial cost component. This gives rise to sudden jump 
in the objective function. As the interstage temperature T^^ 
is increased from a lower value, power for the LP compressor 
increases and the power for the HP compressor decreases. Hence 
a change in the HP compressor corresponds to a downward jump 
while a change in the LP compressor corresponds to an upward 
jump. ; . 

Thus the objective function is discontinuous at 
various points where the size of a compressor changes. At 
these points the function is not differentiable and hence the 
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usual method of optimizations involving derivatives of the 
function cannot be used. For such a case Powell's search 


technique 


(17) 


is used being an ideal choice. 


3.2 OPTIMIZATION TECHNIQUI 


Figure 3.2a shows a flowchart for minimizing a general 

function f(x,/ x,w ..., x_) of n variables: Xw x« , ..., x . 

i- z n 1 2 n 

Considering X = (xj^/ x^ n-dimensional vector, 

the function can also be expressed as f(X). The method uses 
a sequential search technique. In any cycle of minimization, 
a set of n directions D^, D^, . . . , are needed, and to 

start with these are considered as rows of an n x n identity 
matrix. Minimization is first done along each of these n 
directions, as explained later on. Having done that a new 
pattern direction is generated by subtracting the initial 

starting point from the finally reached point X^ (after 
minimizations along each of the n directions). At this point 
the flowchart shows a logical unit to decide whether to retain 
the old set of directions or change them. Thus if either 
f 3 ^ fl or 6 0 (see the flowchart) , the old set of direc- 

tions (D^, D^/ .../ D^) are retained and again minimizations 
are done along each of these directions, unless convergence 
has already been achieved. To the contrary , if both *3 ^ ^ 2 ^ 

and 6 -<1 0 then the old set of directions are changed as 

follows: the direction D along which there had been maximum 

m 

h ''''-It* ' 

decrease in the function value is now replaced by the 

new pattern direction. Having done this now again minimiza- 
tions are done along each of this new set of directions , unless 
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convergence has already been reached. This completes one 
cycle of minimizations , and this cycle is repeated again and 
again until the final points and corresponding function values 
stop differing by a predetermined small quantity. The conver- 
gence is then said to have been achieved and the process is 
terminated, thus getting the minimum of the function f(X). 

Figure 3.2b shows a flowchart for minimizing f(X) 
in a particular direction D, starting from a point This 
problem is equivalent to finding a = a* so as to minimize 
f + cc ^ • Fot any given and ^ one can write, 

f (X a D ) = F(cc), a function of a . 


So the problem is basically a one dimensional minimization of 
F(a) . Quadratic interpolation technique^ ' has been used for 
this purpose. The technique can be well understood by the 
flowchart. In this method a quadratic curve (parabola) is 
fitted through three points of which at least one lies on 
either side of the minimum. Thus if a, b, c be three such 
points (values of a), then a quadratic through these points 
can be taken as ' 


o 

G (^ ) = a^ -t ^2^ -f-a^oc 1 ) 

where the constants a^/ 2 ^ 2 ^ ^3 ^ determined from the 

simultaneous equations 



F(b) 


(3.2b) 



22 


ai + a^c + 330 ^ = F(c) = (3.2c) 

The expressions for a^, a.^/ 83 in terms of a, b, c, 
and Z^ are written in the flowchart (Figure 3.2b). How G(a) 
is considered to be an approxima oion for F(a) and the minimum 
of G(a) lies at 

a* = - 33 /( 283 ) ( 3 . 3 ) 

provided 83 > 0 , which is automatically ensured if at least 
one of a, b, c lies on either side of the minimum point. Now 
cc* is considered as a candidate for the minimum of F(a). If 
the values of F(a*) and G(a*) differ by more than a predeter- 
mined small quantity then one of the three points a, h, c is 
replaced by a* and a new quadratic is fitted through the new 
set of points. The cycle is thus repeated till the conver- 
gence is achieved. 

Flow’chart in Figure 3.2a makes repeated reference to 
flowchart in Figure 3.2b for one dimensional searches. Thus 
the technique essentially consists of a series of one dimen- 
sional minimizations . 

3.2.1 Cost Data: 

For minimizing the objective function for any given 
set of conditions by the technique mentioned in Section 3-2y 
it needs to be evaluated again and again. For evaluating the 
objective function^ costs of compressors^ condensers, elec- 
tricity and water are needed. The results presented in the 
following chapter are based on the costs of compressors and 

^ ^ .... t.. X. 4 .4 X K/f /-I i?* v-' /-I T *i a T Mx:4T/ir "Dpi 1 Vi T 
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and tabulated in Appendix A. Cost of water used is that for 
commercial purposes as obtained from Jal Sansthan, Kanpur and 
it is also given in Appendix A. Cost of electricity for 
various years have been obtained from KESA House, Kanpur and 
are tabulated in Apperdix A. A logistic curve of the form 
a/[b + exp(- cx)] has been fitted through these past data as 
detailed in Appendix B. The choice of such a curve has been 
motivated by the yearly variations in the cost of electricity 
as shown in Figure 3.3. The property values of the refri- 
gerants taken from reference (13) are given in Appendix C. 

3.2.2 Termination Criterion: 

The termination criteria are based on the differences 
between the temperatures and objective function values in 
successive cycles of minimizations. A change in the inter- 
stage or the condensing temperature of 0.01° C or less will 
have no significant effect on any of the quantities of interest. 
Also the annual costs for a 10 ton system are found to be in 
tens of thousands of rupees. So an amount of Rs, 10/- can be 
considered as insignificant in comparison to the total annual 
costs. Thus the termination criteria used during the minimi- 
zation of the cost are as follows: if (i) the temperatures in 
successive cycles of minimizations differ by less than 0.01°C/ 
and (ii) the objective function values in successive cycles 
of minimizations differ by less than 10 (Rs/yr) then the 
process be terminated, thus obtaining the optimum point. The 
corresponding optimal quantities can then be calculated using 
the expressions in Chapter 2, To ensure that the minimum 
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x'eached is the global minimum various starting points were 
tried. In cases involving more than one minim^ the one having 
minimum objective function value was taken as the global minimu 

3.3 TWO STAGE OPTIMIZATION PROCEDURE: 

In the method outlined in the previous section^ the 
minimum of the objective function is reached in a single step. 
However y keeping the discontinuities in the function in view, 
one has to try several starting points selected in the entire 
range. A major computational simplification would result if 
the discontinuities were removed, taking care of them separa- 
tely. This can be achieved by adopting a two-stage technique. 
The first stage involves extensive computer simulation. The 
entire T^^ , T^^ plane is divided into a number of grids and 
performance of the system is studied for each of these grids 
(a set of operating conditions). This study helps determine 
the points where the sizes of the compressors are to be 
changed and the range of temperatures in which each of the 
combination of compressors works. Also, writing down the 
objective function values in the corresponding grids in the 
T , Tj^ plane, a zone of low cost range can be identified as 
explained in Figure 3.4. This zone is the area enclosed by 
the bold lines in that figure. 

The second stage involves finer optimization in the 
grids included in the low cost zone. As the combination of 
compressor sizes which work in each of the grids has already 
been determined in the first step, the cost of these compressors 
can be used for optimizations in the corresponding grids., 
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FIG. 3.4 Detemlnation o£ Low Cost Zones. Numbers in the 
Grids are Costs (Rs/yr) and the Sizes of the Two 
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Thus the costs of compressors remain fixed during optimiza- 
tion in each grid, and the discontinuous element of the objec- 
tive function is removed . 

The costs of the objective function in each of the 
grids are shown in Figure 3.4 together with the combinations 
of LP and HP compressors that would work in each of those 
grids, for a typical case. Three grids have been identified 
there as lying in a low cost zone . But in two of them the same 
combination of compressors works. So finally, only two corioi- 
nations of compressors remain to be tried. Optimization has 
been done with starting points in each of the grids and the 
result has been found to correspond to that obtained by the 
earlier method. The present method, though more time consuming, 
gives a better insight into the problem. 
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CHAPTER - 4 


RESULTS AI4D DISCUSSIONS 


The objective function formulated in Chapter 2 has 
been optimized for various ambient conditions, approach temper- 
atures, degrees of subcooling, degrees of superheat and for 
various tonnage capacities of the system. The results being- 
presented in the following sections pertain to a two-stage 
system for refrigerants R-12 and R-22. Life of the system has 
been taken as 15 years and an operating factor of 0.75 has been 
assumed. In Section 4.5 results are presented for the case 
when the variations in ambient conditions are incorporated as 
discussed in Section 2,2.2. Results of the present study have 
also been compared with available results in the literature. 

4.1 graphical PRESEOTATION of RESULTS: 

Figure 4.1a shows the variation in optimum interstage 
pressure for a 10 ton R-12 system with the evaporator temper- 
ature (T, ) for various inlet water temperatux'es (T .)• Inlet 
water temperatures are kept same as ambient temperatures. In 
this plot an approach temperature of 3°C has been considered 
with no subcooling or superheating. Figure 4.1b shows a 
similar variation for R-22 system. From these figures it is 
possible to read out the value of evaporator 

temperature between -50°C to — 10°C and ambient temperature 
between 15°C to 35°C, for the refrigerants R-12 and R-22.' As 
the evaporator temperature increases, the interstage pressure 
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:j ter oj*timum 3r4aritig between the two stages. Also as 

t--.’ *i‘ * .i*-at temp ’tature is increased mass flow of cooling 

^ 4 1 t^'H-.i tr) irscrease, thus necessitating an upward 

ai tr • -n t in th*r condensing temperature. So for the optimum 

shat in , tc-t xt''‘'*n the two stages P. ^ increases with increased 

^ lopt 

a:n;'i''nt temp/rature. Variations for both R-12 and R-22 are 
sitnilni , th€3 only difference being that the interstage press- 
ur'‘'!t; ar-; greater for R-22 than for R-12 for any given case. 

Figures 4,2a and 4,2b show similar variations in 
OTtimum interstage temperatures for R-12 and R-22 systems, 
res|„>f-fCtively . Whereas in case of interstage pressure the rate 
of increase with increasing evaporator temperature goes on 
increasing^ it shows a reverse trend with the interstage 
aturo. This is due to the relationship between tne 
saturation temperature and saturation pressure for these 
refrigerants . 

Figures 4,3a and 4.3b show the variations in the 
condensing temperature with evaporator temperature having the 
ambient temperature parameter. As the evaporator temper- 
ature is increased the optimum condensing temperature is also 
increased but the increase, is not very pronounced . On the; . 
other hand the condensing temperature is more sensitive to the 
ambient temperature. It is so because the mass flow rate of 
the cooling water is directly dependent on the difference 
between the condensing temperature and the ambient temperature. 
Thus as the ambient temperature increases, the condensing 
temperature is also to. be increased for reducing the cooling 



29 


v.-at-'^r lio-A- the correspording cost. The nature of varia- 

tion tor noth k-12 and R-22 systems are similar. 

Figures 4.4a and 4.4b show the variations in optimum 
cooling wati... r rate with evaporator and ambient temperatures , 
kcepdrip all other parameters fixed. As there is negligible 
variations in optimum operating conditions with respect to 
capacity of the system# the cooling water rate is directly 
proportional to the capacity of the system. For any ambient 
and €:vaporator te^mperatures one can read the value of optimum 
cooling water rate from these graphs. This flow rate should 
be maintained for maintaining the condensing and interstage, 
temperatures at the optimum levels. Hence it is a very 
important quantity. It is evident from the figures that as 
the evaporator temperature is increased the mass flow of 
cooling water decreases. It is so because as the evaporator 
temperature is raised# work of compression reduces so that_ 
there is less heat rejection required. When the ambient temp- 
erature reduces # the water f low required also decreases . It is 
so because in that case the temperature rise of water increases 
so that the same heat is carried by a smaller quantity of 
cooling water. 

Power requirements for LP and HP compressors are 
shown In Figures 4.5a and 4.5b for R- 12 and R-22 systems 
respectively. As the evaporator temperature is raised# the 
power requirement s for both LP and HP compressors decrease. 

The fall is more steep for the HP compressor. The power 
requirement for the HP compressor is seen to be greater than 
that for the LP compressor in all the cases. Further, the 
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difference between the powers required for HP and LP compre- 
ssors is more for lower evaiDorator teraperatures and tends to 
become zero as the evaporator temperature increases beyond -10° 

Values of CO? for optimum operating conditions are 
shown in Figures 4.6a and 4.6'o for refrigerants R-12 and R-22 
respectively. It is seen that for cori'esponding evaporator 
and ambient temperatures, R-12 has a slightly higher COP than 
R-22 when the interstage and condensing temperatures are 
maintained at their optimum levels. 

4.2 EFFECTS OF VARIOUS PARAMETERS: 

4.2.1 Effect of Subcooling and Superheating: 

Effect of subcooling on optimum values of interstage 
pressure, condensing temperature and cooling water rate are 
shown in Figures 4.7a, 4.7b and 4,7c , respectively. It is 
seen that as subcooling of condensate is increased, the optimum 
interstage pressure comes down. Referring to eqn. (2.10b), it 
is seen that as h^ is reduced due to subcooling m^ (the mass 
flow through HP compressor) gets reduced, which means less 
load on HP compressor. To balance it lowered down. 

The optimum condensing temperature increases as the 
subcooling is increased. This is also due to the same reason 
as mentioned above, i.e, the condensing temperature is slightly 
increased to balance the reduced load on the HP compressor. 
Also, as the condensing temperature is increased, temperature 
rise of the cooling water is more and so the mass flow rate 
of cooling water should get reduced with subcooling. This is 
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Table 4.1 shows the effect of superheating of vapour 
in evaporator on optimum interstage and condensing temperature 
for a typical situation. It is seen that with increasing 
degree of superheat both optimum interstage pressure and 
condensir^ temperature increase slightly but the effect is 
negligible . 

For any given subcooling and superheating, the optimum 
quantities can be determined by 

using the following multipliers: 


Pi 


^iopt-^^1- \i- ^ap- ^c- ^s^ 
^iopt^^l' ^wi' ^ap' 


(4.1a) 


T. . -fT, , T , , T , T , T 4- + 273.15 
, _ hopt 1 wi' ap c s^ 

'Th T’* 4T, , T , , t‘ , 0, O-f + 273.15 

hopt 1 wi ap 


(4.1b) 


mw 


m ^4T, , T . , T , T , T ■)• 
wopt 1 wi ap c s 

, T . , T , 0, 04 
wopt 1 wi ap 


(4.1c) 


and 


COP 


COP ^4Tw T T^, TA 

opt 1 wi ap c s 

COP' 4T, , T . , T 0, 04" 

opt^ 1 wx' ap' 


(4, Id) 


The following expressions have been obtained for R-12 
by regression analysis; 


Pi 


Th 


mw 


1 - 0.01654 - 0.00661 + 0.000105 T^ (4.2) 
1 + 0.000256 T + 0.000175 - 0.0000125 (4.3) 
1 - 0.01701 T , + 0*003257 - 0.000210 tJ (4.4) 


j 3 
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Analogous expressions tor the case of R-22 system are: 


F . 
Pi 

= 1 - 0,1010 T - 0.00377 tJ + 0.000104 T^ 

G G C 

(4.6) 

T 

Th 

= 1 + 0.000188 T + 0.000106 T^ - 0.0000117 T^ 

G C G 

(4.7) 

F 

mw 

= 1 - 0.01569 T^ 4- 0.00251 T^ - 0.000157 T^ 

c c c 

(4.8) 

F 

COP 

= 14- 0.01001 T^ - 0.00211 T^ 4- 0.000169 T^ 

c c c 

(4.9) 

The 

effect of superheating has been neglected. Thus, 

first 


the optimum quantities are to be found from the graphs given 
in Figures 4.1 to 4.6 for the case without subcooling and 
superheating , and then these quantities are to be multiplied 
by the corresponding multipliers given above, 

4,2.2 Effect of Approach Temperatures: 

In the above sections all the results nave been 
presented for 3°C approach. But/ in individual cases they 
may be required to be different than 3°C. So the effect of 
various approach temperatures on the optimum quantities have 
been considered. It has been found that the optimum values 
of interstage pressure, condensing temperature and cooling 
water rate all increase as the approach temperature is 
increased . But optimum value of cop reduces with increasing 
approach temperatures, fiowcver, the effect of approach temp- 
erature is more pronounced on cooling water rate than the 
other quant.ities. It is so because the temperature rise of 
the cooling water id directly related to the approach temper- 

Figure 4.8 quantitatively shows the effect of approach 
on the cooling water rate. 
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For apnroach ouv-*r than 3°C, multipliers defined by: 


■pi 


P 

^iopt^‘l' 

p.' 4 t, 


T . , 

wx 


T 


ap^ 


r ^ 
c 


T 4* 
s 


« . \ * "1 

lopt 1 wl 


3 / T . 


T > 
s 


( 4.10a) 


Th 


nopt 1 V.X ap c s 


"^noot ' ' 1 ' * vv’i 


fT, , T_. , 3, T^, + 273.15 


mw 


"'ap- ‘s» 


COP 


COP .,4T, , T , , T , T , T .)■ 

opt 1 wi ap c s 

COP' ^-fT, , T . , 3, T , T ■)• 
opt 1 wd G s 


(4.10b) 


(4.10c) 


(4.10d) 


have been determined by regression analysis. For R- 12 system, 
these are: 


Si 

= 

0.96830 - 0.0001363 T 

ap 

+ 0.002106 

ap 

(4.11) 

°Th 

s 

0.95736 - 0.0005398 

ap 

+ 0.002380 tI 
ap 

(4.12) 

G 

mw 

= 

0.80860 + 0.09529 - 

ap 

0.00938 

(4.13) 

Sop 


1.02809 - 0.001064 

op 

- 0.001938 

ap 

(4.14) 

and 

for.. 

R-22 system. 



Si 

= . 

0 .9914 4 - 0,002 563 4 

ap 

+ 0.001407 

(4.15) 

’^Th 


0.98296 - 0.0033177 

ap. 

+ 0.002476 

(4.16) 

G 

mw 


0.78609 + 0.10942 T ^ - 

ap 

0.012483 

ap 

(4.17) 

Sop 


1.01943 + 0.00340 - 

ap 

- 0.002300 

ap 

(4.18) 
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Thus to compute the various quantities tor approach 
temperature dirtf^rent than 3°C, one has to multiply the results 
with 3°C approacn by the corresponding multipliers given above. 

4.3 ILLUSTRATIOf.': 


It is to select the optimum 2-stage R-12 vapour-compre- 
ssion refrigeration system for a given set of requirements and 
to recommend the optimum operating conditions. 

Available Information: 

1. Capacity of the plant: 10 tons 

2. Evaporator temperature required; -35®C 

3. Ambient temperature; 30° C 

4. Approach between the condensing and exit cooling water 

temperatures ; 5° C 

5. Degree of subcooling of the condensate possible; 5°C. 


Solution; 

First the results are found for the case without sub- 
cooling and superheating and with approach temperature of 3°C. 
From Figure 4.1a, 


Optimum interstage pressure, 


3.5 bar 


From Figure 4.2a, 


T’ = 4.2°C 

lopt 


From Figure 4,3a, Optimum condensing temperature, = 48.5 ( 
Prom Figure 4.4a, Optimum cooling water rate, " 0.853 Kg/s 
Prom Figure 4.5a, 


Power required by LP -compressor 


12.3 HP, and 
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Power required uy liP -compressor = 15.5 HP. 

So both LP esnd '’-IP comrressors need to be of 15 HP capacity^ 
since the capacity next higher to 10 HP is 15 HP. 


From Figure 4.6a/ 

COP of the system for optinmra conditions = 1.73 5 


To correct for approach temperature, subcooling and super- 
heating the multipliers given by eqns. (4. 2-4. 5 and 4.11-4.14) 
ax'e computed as : 


F„. = 1 - 0. 01554 X 5 - 0.00661 x 5^ + 0.000105 x 5^ 

Pi 

= 0.76517 

= 1.00409, = 0.97012, = 1.032375 

Gp. = 0.96830 - 0.0001363 x 5 + 0.002106 x 5^ = 1.02027 

°Th = = 1.05055, G^^jp = 0.97432 

So the corrected values of various quantities are; 


P- ....L. 
lOpt 

i .e . / 


P^ ^ 

lopt 


X F 


Pi 


X G 


Pi 


3,5 X 0.76517 x 1.02027 


P. ^ 
lopt 

T 

hopt 

(Note; 


= 2 ,630 bar 

= '^hopt + Sh - 

The multiplies for the condensing temperature are 
obtained in terms of absolute temperatures.) 


So, 


T 

hopt 


54.33«C, 


ih _ = 0.86934 Kg/s 

wopt 


COP 


opt 


1.745 


"Phus , one selects ; 
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X ,5 ^ ™ 2«630 osjt 
Gorki*; ns in,} t-nin.'^^'racure = 5-4-38®C 

Cooling wat^r rata = 0. 06^*34 Kg/s 

capacity of Ln-couK;reasor = i5 r.r 
.capacity or -compressor ■ 


15 HP. 


4.4 cnHPAHl30i4 VriTH AVAILABLE DATA; 

in the present study, both the condensing and the 
interstage temperature have been considered as decision 
variables and their optimum values have been found together 
to minimize the total annual cost of the system. But th 
results available in literature have all been obtained for 
fixed condensing temperature with only interstage pressure as 
the decision variable. Thus to ensure a meaningful comparison., 
the condensing temperature has been fixed and one aimensional 
minimization has been carried out with respect to interstage 
temperature. The results, thus obtained, are shown in Table 
4.2, together with results available in the literature. It is 
seen that the optimum interstage pressure is about 9% m 

,, . „.,-c « ■ »■>•= •• 

the geometric mean pressure. The difference betwe 
present study and other investigators is seen to be small 
for = 42.9»C. But, for = 55.5”C, the difference is 
significant. It is 694 more as compared to the result of 
verma et and 4.W more as compared to that of Prasad 

et al.^^®^ 

The optimum interstage pressure by the present method 
= seen to be greater than that of all other investigators. 
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It is i-.ft-dnly due to the s.-dtt in the optimization criterion. 

4.5 i ATICdd id TFAP^RATuEE : 

Thv- eiifeot or year-round variations in the ambient 
temperature haa oec.n incorporated as discussed in Section 2.2. 
Figures 4.9a, 4,ijb, 4,9c and 4.9d show respectively the varia- 
tions of optimum values of conciensing temperature, interstage 
pressure, powe?r requirements of LP and HP compressors and 
cooling water rate, with the evaporator temperature. The 
nature of their variations is similar to those for the case 
when the arrbient temrierature is considered to be a fixed desig, 
value. There is a difference between the two cases only in 
numerical values as can be seen by comparing the graphs for the 
two cases. Table 4,3 shows a comparison between the two cases 
for a 10 ton 2-3tage R--12 system. The maximum ambient temper- 
ature has been considered as the design ambient temperature 
for the first case. It is seen that when the effect of the 
variations in ambient temperature is considered , both the 
cooling water rate and the annual cost comes down. Thus, for 
the evaporator temperature of -30°C, -40°C and -50°C, savings 
of Rs.1982/-, Rs. 6928/- and Rs.8917/- result every year. 
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CHAPTER - 5 

CONCLUSIONS AND SUGGESTIONS 


5.1 CONCLUSIONS; 

From the present study the following conclusions are 
arrived at: 

1. There is a significant difference between optimal systems 
obtained from thermodynamic and cost considerations. The 
total cost of the system comprising initial and running 
expenditures for the entire life (in terms of present 
value) has been considered as the criterion of optimi- 
zation . 

2. The optimum interstage pressures found in the present 
study are greater than those found by all other investi- 
gators. The present value is found to be about 15% 
higher than the geometric mean pressure. 

3. The condensing temperature has also been considered as a 
decision variable because the cost of cooling water is 
directly related to it. An optimum condensing temper- 
ature , together with economic cooling water rate has 
been found out. 

4. Multipliers have been determined to account for the 
effects of subcoolingy superheating and various approach- 
temperatures ^ within insignificant deviations. The 
optimum interstage pressure decreases with subcooling 
and slightly increases with superheating in the evapo- 
rator. On the other hand, the optimum condensing 
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temperature increases ^ile the economic cooling water 
rate reduces as the subcooling is increased . The varia- 
tion in approach temperature has maximum erfect on 
economic cooling water rate. 

5. The power requ.irement for HP -compressor has been found to 
be more than that for the LP-compressor . The difference 
between these power requirements is more for lower evapo- 
rator temperatures and tends to zero as the evaporator 
temperature is increased beyond -10° C. 

6. The system has been optimized both for a fixed design 
ambient temperature and incorporating the effect of year- 
round variations in the ambient temperature. The system 
optimized for the latter case has lesser total cost as 
compared to that for the first case. Thus, even though 
the system is not optimal for any particular ambient 
temperature, it is optimal for the year as a whole. 

7. The expressions and methodology developed are general 
and can be used for any set of ambient conditions and 
prevailing costs. Also the method can be used for any 
refrigerant for which the correlations for the properties 
are available . 

5.2 SUGGESTIONS: 

The best system would be the one which is optimal for 
all ambient temperatures round the year. But the 
components of the system has to be selected to meet the 
the worst airbient temperature. However, as the ambient 
temperature varies, the condensing temperature and the 


1 . 
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interstage teraparature v.’ill change. So a control device. ^ 
would be required to raaintain the optimurri interstage 
tempx'jrature for each ambient temperature, 

2. In tiie present study trie cooling load has been considered 
to be fixed. But as the ambient temperature varies the 
load also varies. The present work can be extended to 
include the effect of variations in the cooling load with 


time . 



41 


REFERENCES 


(1) Threlkeld , J.L., "Thermal Environmental Engineering"/ 
Prentice Hall Inc., Englewood Cliffs, N.J,, USA, 1962, 

P . 67 . 

(2) Soumerai, Henri P., "Simplified Analysis of Two-stage 
Low Temperature Refrigeration Systems" , Refrigerating 
Engineering, Vol. 61, 1951, pp. 746-50 & 802. 

(3) Keshwani, H.B. and Rastogi, K.V., "Optimum Intercooler 
Pressure in Multistage Compression of Refrigerants" , 

All India Symposium on Refrigeration, Air-Conditioning 
and Environmental Control in Cold Storage Industry, 
I.I.T., Kanpur (India), 1968, pp. E6.1-13. 

(4) Lai, Ram, "Selection of Intermediate Pressure in 
Compound Refrigeration System", M.E. Thesis, Dept, of 
Mech. Engg., Univ. of Roorkee (India), 1970. 

(5) Verma, H.K. and Charan, V., "Intermediate Pressure for 
Two-stage Compression System with Flash Intercooling" , 
Second National Symposium on Refrigeration and Air- 
Conditioning, Univ. of Roorkee (India), 1973, pp. 75-79. 

(6) Gosney, W.B., "Compound Compression Refrigeration 
Cycles", The Institute of Refrigeration (London), 1966-67, 
pp. 1-11. 

(7) Arora, C.P. and Dhar, P.L., "Optimization of Multi-stage 
Refrigerant Compressor" , Progress in Refrigeration Science 
and Technology, XIII International Congress of Refriger- 
ation (HR), 1973, pp. 693-700. 

(8) Prasad, M., "Optimum Interstage Pressure for Two-stage 
Refrigeration System", ASHRAE Journal, Vol. 23, No. 1, 
1981, pp. 52-60. 

(9) Gupta, v.K. and Prasad, M., "Correlations for Optimum 
Interstage Pressures for R-22 Two-stage Refrigerating 
System Having Effects of Subcooling and Superheating" , 
Proceedings of 4th ISME Conf . on Mech. Engg., Mech. 

Engg. Dept., Univ. of Roorkee (India), Oct. 19-20, 1981, 
pp. 495-499. 

(10) Agrawal, H.C. and Kapoor, A.K., "Design and Development 

Solar water Heaters for Indian Conditions (Phase-1, 

A Preliminary Study)" , April, 1975, Table 2.3, pp. 2.5- 

2 . 8 . : 

(11) "Refrigeration and Air-Conditioning Data Book, Design 
Volume", 1957, p. 18.02. 



42 


(12) Prasad, S., "Optimum Inter-stage Pressure for R-12 
Refrigeration System", Institution of Engineers (India), 
Vol. 12, Part S4, ST March, 1979, pp. 118-120. 

(13) Prasad, M., "A. Note on Refrigeration and Air-Condi- 
tioning", Dept, of Mech. Engg. , I .1 .T, , Kanpur , pp. 
4.67-4.69 \i97B9, 

(14) Borbely , P., "Economic Insulation Thickness", XV Inter- 
national Congress of Refrigeration, Venezia, 23-29 Sept. 
1979, Dl-32. 

(15) Prasad, M. and Gupta, V.K., "Economic Cooling Mater 
Rate for Two-stage R-12, R-22 and R-717 Systems", 
Proceedings of Seventh Biannual National Symposium on 
Refrigeration and Air-Conditioning, I.I.T. Kharagpur, 
December 22-30, 1980, pp. 15-20. 

(16) Prasad, M. and Kaul, P.N., "A Short Term Course on 
Advances in Refrigeration and Air-Conditioning", 
Organised by Dept, of Mech. Engg., I.I.T. , Kanpur, 

May 11-31, 1981, pp. 1.39-1,40. 

(17) Powel, M.J.D., "An Efficient Method for Finding the 
Minimum of a Function of Several Variables without 
Calculating Derivatives", Computer Journal, 7(4), 1964, 
pp. 303-307. 

(18) Fox, R.L., "Optimization Methods for Engineering Design" 
Addison Mesley, Reading, Mass, 1971. 

(19) ■ Prasad, M. and Prasad, S., "Optimum Interstage Pressure 

for R-12 System Including Effects of Subcooling and 
Superheating" , Journal of the Institution of Engineers 
(India), Vol. 63, Nov. 1982, pp. 112-116. 



APPENDIX - A 


A.l Cost o£ Compressors for Refrigerant R-12: 


1 

Capacity ; 
(HP) ; 

1 r 1 t 

0.5 : 1.0 • 2.0 i 3.0 ; 

t I t i 

' 

5.0 

7.5 i 
1 
1 

10.0 

i I 

; 15.0 ; 20 

1 1 

2050 3260 4850 

8270 

17800 

26400 

55900 75SC 

A. 2 Cost 

of Compressors for R-22 

; 




Capacity 

1 1 1 

1 « 1 

; 3 ; 5 1 7 . 5 

t 1 

I t 

: 10 : 

15 

t 

; 20 

(HP) 

t 1 I 

1 

1 

t 


1 

1 - 

Cost 

4800 6500 7500 

18000 25000 

4 5000 

(Rs. ) 






A. 3 Cost 

of Mater for Kanpur City: 




Year 

1 

f 

1 For House Use 

1 

1 

» 

1 

1 

1 

1 

1 

For Industrial Use 

1970-1978 

3 

Rs . 0.2 5/m 


Rs 

. 0.35/ra^ 

1979-1982 

3 

Rs, 0,65/rn. 


RS 

. 1.65/m^ 

A. 4 Cost 

of Electricity: 





Year 

1 1972 i 1976 i 1977 

1978 

i 1979 

1 

t 1 . , 

i 1981 i 1983 


0.36 0.43 0.45 0.50 0.55 0.60 0.65 


Cost 

(RsAWh) 
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APPETOIX - C 

PROPERTIES OF REFRIGERANTS 

C.l Refrigerant R-12: 

c_(T) = 0.59524 + 0.0018175 (T + 15) 

P 

s^(T) = 0.142093 + 0.33747 T/lOO - 0.0399606(T/100)^ 

+ 0.019869(T/100)^ + 0.014023 (T/lOO) ^ 

h (T) = 188.86 + 0.440278 T - 7 .02007(T/100)^ - 

y 

5.0765l(T/100)^ - 3.82545(T/100)^ 

h^CT) = 36.1554 + 0.928108 T + 6 . 89916 (T/lOO) ^ 

+ 3.73414(T/100)^ + 5.91673 (T/lOO) ^ 

p(T) = 3.08968 + 0.10162 T + 12 .6413(T/lOO)^ 

+ 6.74792(T/100)^ + 1 .01057(T/100) ^ 

C«2 Refrigerant R-22: 

2 

c (T) = 0.60 + 0.0005246(T + 45 + 0.043686(T + 45) ) 

p 

• for -60 T 10 

= 0.70114 + 0.0029529(T - 10 + 0.010472(T - 10)^) 

for 10 T 50 

s^(T) = 0.181132+0.43667 T/lOO - 0.0596756(T/100)^ 

- 0.0234343(T/100)^ +0.0638047 (T/100)^ 

h (T) = 251.106 +35-4577(T/iOO) - 19 . 4993 (T/lOO) ^ . 

- 7.62005(T/100)^ - 11 . 67 56(T/100) 
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h-(T) = 46.2102 + 1.20394 T + 6.828(T/100) 

- 8.99387(T/100)^ + 20.06l2(T/100) 


p(T) 


4.98105 + 0.151746 T + 19 . 5655(T/100) 
+ 9.8434(T/100)^ + 2 .22194(T/100) . 



TABLE 4.1 

Effect of 

Superheating on 

and 

T 

hopt 

°c i 

P . ^ , bar 

i 

' T . °C 

1 

1 


s J 

lopt 

1 hopt' 

( 

1 


0 

3.303 

53.07 

Refrigerant: R-12 

5 

3.308 

53.08 

Wl 

= 35°C 

10 

3.310 

53.08 

^c 

= 0°C 

15 

3.322 

53.09 

T 

ap 

= 3°C 

20 

3.328 

53.10 

^1 

= - 40° C 


TABLE 4.2 Comparison of Optimum Interstage Pressures from 
Different Sources 


Operating 

Conditions 

(“C) 

T ' T 

^1 i ^h 

Present 

Method 

bar 

1 

Ref .(19) 

bar 

Ref . (5) 

bar 

Ref . (7.) 

bar 

- - J 

Ref . (12) 

bar 

1 , - - - - 

(Pi,P, 

bar 

-11.7. 

42.9 

5.043 

5.02 

5.10 


5.06 

4.62 

-11.7, 

55.5 

6.112 

5.85 

_ 

5.66 

5.99 

5.34 


TABLE 4.3 Comparison Between Results for Design Ambient Ternpe: 

ature and with Variations in Ambient Temperature 


Operating 

variables 

With Effect of Varia-lWith Design 
tions in Ambient ; ^ 

Temperature I wmax 

* ... - 

Conditic 

= 40° C 


t 1 

I'm''' \ 

1 : 

V- ■ 

T , 
ap' 

m , / ; 

wopt , 

^opt ' 

i ^ i 

I wopt I 

^opt' 

tons 

t ,* 

t t 

°c 

°C 

Kg/S 1 

-I 

Rs/yr 

: Kg/s 

Rs/yr 

10 

-30 

0 

3 

0.7901 

85323 

0.8821 

87305 

10 

-40 

0 

; 3; V 

0.8928 

101713 

0.9835 

108641 

10 

-50 

0 


1.0492 

123242. 

1.2053 

132159 


r — I J ^ 

r : 


ii mi 


(«n+a) 


/' T~^f 

f r> 

1 jf 


(-W 5 I 


^ / lil^^UA: 

§m- 2 } 4 b 

'^nr-^iiii 


) ^ 


FIG-S‘l SCHEMATIC DIAGRAM OF 
71 -stage SySTEM 
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FIG‘ 3-3 f>-/> DIAGRAM FOR 71' STAGE SYSTEM 
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fig b-6 a real condensation process 
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